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Summary 


This  lecture  is  devoted  to  the  basic  physical  properties  and  the  calculation  methods  of  the  artificial 
cavitation  (ventilation)  flows.  Main  attention  is  paid  to  the  three-dimensional  ventilated  supercavities  past 
the  disk  and  blunted  cavitators. 

Derivation  of  approximate  equations  of  the  axially  symmetric  supercavities  with  using  the  theorem  of 
momentum  and  also  the  asymptotic  theory  of  slender  body  is  presented. 

The  G.V.Logvinovich  independence  principle  of  the  cavity  expansion  is  stated  and  grounded.  It 
gives  a simple  and  general  method  of  investigating  the  unsteady  problems. 

Calculation  of  the  ventilated  cavities  has  some  peculiarities  caused  by  gravity  effect,  presence  of  the 
cavitator  angle  of  attack  and  taking  account  of  balance  between  the  gas-supply  into  the  cavity  and  the 
gas-leakage  from  the  cavity.  They  are  discussed  here  as  well. 


Introduction 


The  basic  similarity  parameter  of  cavitation  flows  is  the  cavitation  number  [1-4]: 


o = 2[p--fc) 

pv: 


(i) 


where  px  is  the  pressure  at  infinity;  p,  is  the  cavity  pressure;  p is  the  water  density;  Vx  is  the 
mainstream  velocity.  When  a body  moves  horizontally  on  the  immersion  depth  H,  we  have  px  = pgH  , 
where  g is  the  gravity  acceleration.  The  supercavitation  flow  regime  corresponds  to  small  magnitudes 
a < 0. 1 . When  the  magnitude  of  a decreases,  the  supercavity  dimensions  increase. 

When  p = const,  there  are  three  possibilities  of  the  a reduction: 


1 ) Increasing  the  mainstream  velocity  Vx . With  decreasing  the  mainstream  velocity  the  rarefaction 
past  bodies  increases,  and  cavities  are  created  by  natural  way  when  the  pressure  decreases  to  a value 
close  to  the  saturated  water  vapor  pressure  pv=  2350  Pascal  (at  the  temperature  20 °C).  Such 
supercavities  are  called  the  natural  or  vapor  supercavities.  The  parameter 


c. 


2£P~  ~ Pr ) 


(2) 


is  called  the  natural  or  vapor  cavitation  number.  The  natural  supercavities  are  filled  by  water  vapor  under 
pressure  pc=pv.  The  natural  cavitation  regime,  when  av  < 0.1,  is  attained  when  the  mainstream 
velocity  is  very  high:  V^>  50  m/sec. 


Paper  presented  at  the  RTOAVT  Lecture  Series  on  ", Supercavitating  Flows  ”,  held  at  the  von  Karman 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 


11-2 


2)  Decreasing  the  difference  Ap  = - pc  owing  to  decreasing  the  ambient  pressure  p„ , for 

example,  in  closed  cavitation  hydro-tunnels  [4,  5],  The  super-cavitation  flows  in  the  stream  with 
artificially  reduced  pressure  are  sometimes  called  the  artificial  ones,  although  they  do  not  distinguish  on 
natural  vapor  supercavities  by  their  physical  properties. 

3)  Decreasing  the  difference  Ap  = px-  pc  owing  to  increasing  the  cavity  pressure  pc . This  is  easy 
attained  by  blowing  the  air  or  another  gas  into  the  rarefaction  zone  past  a non-streamlined  body.  Such 
supercavitation  regime  is  called  the  artificial  cavitation  or  ventilation.  The  method  of  obtaining  the 
ventilated  cavities  was  proposed  independently  by  L.A.Epshtein  and  H.Reichardt  in  1944  - 1945. 

The  ventilated  cavitation  regime  o <0.1  is  easy  obtained  at  the  significantly  lower  mainstream 
velocities  ~ 10  m/sec.  That  is  why,  it  is  widely  applied  in  bench  investigations  of  the  supercavitation 
flows.  The  Fig.  I presents  a typical  scheme  of  an  experimental  apparatus  to  investigate  the  ventilated 
cavities. 


Fig.  1 . Scheme  of  apparatus  to  create  a ventilated  supercavity 


1 Main  problems  of  artificial  cavitation 


It  follows  from  the  similarity  theory  of  hydrodynamic  flows  [3]  that  the  shape  and  dimensions  of  both  the 
natural  vapor  supercavity  and  the  artificial  ventilated  cavity  must  be  the  same  when  the  cavity  number  is 
equal.  Flowever,  in  reality  the  full  similarity  does  not  fulfil  due  to  a number  of  causes. 

The  main  peculiarity  of  the  artificial  cavitation  flows  is  a consequence  of  comparatively  small  values 
of  the  velocity  V=  10-4-  100  m/sec.  In  this  case  the  Froude  number 


Fr  = 


(3) 


where  Dn  is  the  cavitator  diameter,  has  moderate  magnitude. 

This  means  that  the  gravity  effect  is  considerably  more  for  the  artificial  cavity  when  the  cavitation 
number  is  the  same.  It  causes  to  the  more  essential  deformation  of  the  artificial  cavity  shape.  The 
floating-up  the  ventilated  cavity  tail  may  be  essential  so  that  it  completely  determines  the  cavity 
closure  character  and  type  of  the  gas  leakage  from  the  cavity. 

Thus,  the  Froude  number  (3)  together  with  the  cavitation  number  (1)  is  the  basic  similarity  parameter 
for  the  artificial  supercavitation  flows. 
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1.1  Methods  of  creation  of  the  ventilated  super-cavities 


We  know  the  following  methods  of  creating  the  artificial  ventilated  cavities  (Fig.  2): 

1)  Method  of  a gas  source  [6]  consists  in  creating  a gas  envelope  around  a body  by  the  air  jet  ejection 
from  the  body  nose  opposite  to  the  stream  (Fig.  2,  a).  It  is  obvious  that  in  this  case  the  gas  pressure  must 
exceed  the  pressure  of  the  water  in  the  nose  stagnation  point.  From  the  Bernoulli  equation  [7],  we  obtain 
estimation  of  the  necessary  velocity  of  the  air  Vg  : 


V > V 

g 00 


-£-(1+<t)®28.6 

\pg 


f 

V 


(4) 


where  p is  the  air  density.  One  can  see  that  the  gas  jet  velocity  and,  hence,  the  gas  rate  must  be  very 
high. 

Owing  to  the  high  velocity  of  the  gas  flow  in  the  cavity,  the  cavity  walls  are  strongly  perturbed  (the 
cavity  is  nontransparent).  The  flow  represents  a strongly  turbulized  two-phase  zone  at  the  tail  cavity  part 
and  in  the  wake. 

This  way  of  the  artifitial  supercavity  creation  is  not  applied  in  practice  due  to  its  lacks. 

2)  L.l.Sedov  scheme  [7]  consists  in  ejecting  the  slender  water  jet  from  the  body  nose  opposite  to 
the  stream  with  simultaneous  air-supply  into  the  stagnation  zone  (so  called  "jet  cavitator",  Fig.  2,  b).  This 
flow  scheme  represents  the  inverted  cavitation  flow  by  Efros-Gilbarg  with  the  reentrant  jet  [7],  The  nose 
stagnation  point  is  displaced  into  the  stream,  and  in  the  ideal  case  of  complete  recovery  of  the  pressure  in 
the  tail  part,  the  body  must  suffer  thrust 

T = pS/ViV„,  (5) 


where  S - is  the  jet  section  area;  V j is  the  average  velocity  in  the  jet. 

In  real  flows,  it  is  impossible  to  attain  smooth  closure  of  the  cavity  on  the  body  tail  and  complete 
recovery  of  the  pressure.  Therefore,  the  arising  thrust  will  be  smaller.  Nevertheless,  this  scheme  is  very 
perspective.  When  <7  — > 0.  we  obtain  the  infinite  Kirchhoff s cavity.  In  this  case  the  body  with  the  jet 
cavitator  suffers  half  drag  in  comparison  with  the  body  with  the  solid  cavitator. 


Fig.  2.  Ways  of  creation  of  ventilated  cavities:  a - gas  source;  b - L.I.  Sedov's  scheme; 

c - blowing  past  a sharp  edge 
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3)  Gas-supply  past  a sharp  edge.  The  most  stable  and  smooth  surface  of  the  cavity  is  observed,  if  the 
cavity  separation  occurs  on  the  sharp  cavitator  edge  with  simultaneous  gas-supply  to  the  separation  place 
(Fig.  2,  c).  it  is  possible  to  supply  gas  in  any  point  of  the  body  after  the  fully  development  of  the 
supercavity  (so  called  "distributed  gas-supply").  Experiments  showed  that  location  of  the  place  of  the 
gas-supply  into  the  steady  cavity  has  no  essential  significance. 

The  way  of  creation  of  ventilated  cavities  past  cavitators  with  the  sharp  edge  (for  example  a disk)  is 
the  most  frequently  used. 


1.2  Coefficient  of  gas-supply  rate  into  the  cavity 


The  main  problem  for  the  ventilated  cavitation  flows  is  a problem  on  the  gas  quantity  necessary  for 
supply  into  the  cavity  to  maintain  the  supercavity  of  given  dimension  (i.e.  to  attain  the  given  cr ).  In  the 
case  of  axisymmetric  cavitator,  the  dimensionless  coefficient  of  the  rate  is  introduced  [2,  3]: 


<2  = 


Q 


(6) 


where  Q is  the  volumetric  gas  rate  at  the  pressure  pc , its  dimensionality  is  m 3 / sec. 


Obviously,  the  gas-supply  rate  Qin  must  be  equal  to  the  gas-leakage  rate  Qout  when  motion  is  steady. 
For  unsteady  ventilated  cavities,  the  gas  leakage  from  the  cavity  is  the  unknown  time  function  depending 
on  the  variable  pressure  in  the  cavity  pc  (?)  and  on  the  gas-leakage  type.  Determination  of  the  function 

Qout  (?)  for  concrete  conditions  represents  a basic  problem  of  the  artificial  cavitation  theory. 

A unified  theory  of  the  gas  leakage  from  the  ventilated  cavity,  that  allows  to  calculate  the  function 
Q((J)  for  any  values  of  the  flow  parameters,  does  not  exist.  The  mechanism  of  the  gas-leakage  may  be 
considerably  distinguished  in  dependence  on  the  parameters  Fr,  <7,  <r„.  The  gas-leakage  rate  is  also 
different  in  the  cases  of  free  cavity  closure  and  the  cavity  closure  on  a body.  Approximate  estimations 
were  obtained  for  some  types  of  the  gas-leakage. 


1.3  Hydrodynamic  schemes  of  gas  leakage  from  the  cavity 


Observations  show  that  two  principally  different  mechanisms  of  the  gas  leakage  from  the  stable  artificial 
cavities  exist  when  the  cavity  closure  is  free  [2,  3], 


Fig.  3.  Scheme  of  the  portion  gas-leakage 

When  magnitudes  of  the  numbers  Fr  and  a are  high,  i.e.  the  gravity  effect  is  not  considerable  and 
the  cavity  shape  is  close  to  the  axisymmetric  one,  the  cavity  tail  is  filled  by  a foam  which  is 
periodically  rejected  by  portions  in  the  fonn  of  toroidal  vortices  (so  called  the  portion  leakage  or  the 
first  type  of  gas  leakage  from  the  cavity.  Fig.  3). 

When  Fr  is  fixed  and  <7  is  small,  i.e.  the  gravity  influence  is  considerable,  the  flow  in  the  cavity 
end  becomes  well-ordered.  The  cavity  ends  by  two  hollow  vortex  tubes  carrying  over  the  gas  from  the 
cavity  (so  called  leakage  by  vortex  tubes  or  the  second  type  of  the  gas  leakage  from  the  cavity.  Fig.  4). 
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Besides  the  first  and  second  types  of  the  gas  leakage  from  stable  cavities,  the  third  type  of  the  gas 
leakage  exists  owing  to  pulsation  of  the  unsteady  ventilated  cavities  [8,  9], 


C 20^0 
OOQ 


Fig.  4.  Scheme  of  the  gas-leakage  by  vortex  tubes:  a - view  from  side,  h - view  from  above 

1.3.1  The  first  type  (portion  gas-leakage) 


The  first  type  of  the  gas  leakage  from  the  artificial  cavities  is  not  well  studied.  This  is  explained  by  that  in 
this  case  the  processes  at  the  cavity  end  are  considerably  non-stationary  and  multi-parametric  ones.  It 
was  shown  [3,  10]  that  they  are  determined  by  not  only  the  similarity  parameters  o,  Fr  but  depend  on 
both  the  Reynolds  number  Re  and  the  Weber  number  We: 


Re 


We 


a 

2 pVlDn  ’ 


(7) 


where  v is  the  kinematic  coefficient  of  viscosity;  a is  the  coefficient  of  water  surface  tension.  They  also 
depend  on  degree  of  the  free  stream  turbulence  and  other  external  perturbations: 

Q=F(a,av,Fr,Re,We,...).  (8) 


It  was  shown  experimentally  that  decrease  of  numbers  Re  and  We,  i.e.  increase  of  viscosity  and  surface 
tension  of  liquid,  causes  to  decrease  of  the  gas  leakage  rate  from  the  cavity.  Some  data  and 
hypothesises  about  the  first  type  of  the  gas  leakage  are  presented  in  the  papers  by  L.A.  Epshtein  [10, 
11], 

Qualitatively,  in  this  case  processes  may  be  described  in  the  following  way.  Fluid  particles  flowing 
around  the  axisymmetric  cavitator  in  the  radial  planes  collide  in  the  cavity  closure  zone  and  create  a 
non-stationary  unstable  reentrant  jet.  When  the  reentrant  jet  is  distorted,  the  two-phase  medium  (foam) 
is  formed.  It  partially  fills  the  cavity  and  moves  within  the  cavity  in  the  form  of  a toroidal  vortex.  The 
outer  part  of  this  vortex  moves  along  the  flow  direction  under  action  of  viscous  shears  on  the  cavity 
boundary.  Its  central  part  moves  opposite  to  the  stream  under  action  of  the  pressure  gradient.  The 
vortex  continuously  grows  owing  to  new  foam  portions.  The  foam  rejection  occurs  in  the  form  of  a 
ring  toroidal  vortex  or  its  part,  or  separate  foam  portions  (see  Fig.  3),  when  the  friction  forces  become 
equal  to  the  reentrant  motion  momentum.  Also,  the  continuous  foam  leakage  from  the  zone  past  the  rear 
stagnation  point  exists  simultaneously  with  the  portion  leakage. 

In  the  work  [10]  was  shown  using  the  methods  of  the  similarity  theory  that  the  expression  for  a 
dimensionless  period  of  the  foam  portions  rejection  should  be  in  the  form: 


T = 


D„ 


(9) 


where  cx  is  the  cavitation  drag  coefficient,  g is  the  friction  coefficient,  BX,B2  are  the  function  of  Re, 
We. 
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It  follows  from  the  experimental  data  [12]  that  an  evaluation  of  Strouhal  number  Sh  for  this  process 
is  : 

Sh  = ^ :-«0.315, 

V- 

where  Lm  is  the  average  cavity  length.  The  real  process  frequency  / has  values  10  4-100  Hz  for  various 
tests.  It  is  possible  to  observe  the  portion  gas-leakage  with  a stroboscope  or  using  photography  with  the 
short  exposure  time  10-4-  10~5  sec. 

In  the  case  of  the  flow  close  to  the  natural  vapor  cavitation  when  the  gravity  effect  is 
insignificant,  G.V.Logvinovich  established  the  semi-empirical  law  of  the  gas-leakage  [2]: 

Q=rv~sc\y-A  (io) 

where  y = 0.01-^0.02  is  the  empirical  constant,  Sc  is  the  area  of  the  cavity  mid-section. 


0.02  0.06  a 

Fig.  5.  Influence  of  the  Froude  number  on  Q (cr) : 
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Fig.  6.  Experimental  dependence  Q(<7) 
Fr=  13.3 


1 - Fr  = 19.3;  2-16.5;  3-14.6;  4-12.7;  5-11.0 
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The  portion  type  of  the  gas-leakage  with  the  periodically  arising  and  distorting  reentrant  jet  is 
characteristic  also  in  the  case  of  supercavitation  flow  around  hydrofoils  with  great  aspect  ratio  [6]. 


1.3.2  The  second  type  (gas-leakage  by  vortex  tubes) 


When  the  second  type  of  the  gas  leakage  from  the  artificial  cavity,  the  process  is  stationary  and 
completely  determined  by  the  gravity  effect  on  the  cavity.  The  second  leakage  type  is  the  most  studied 
both  theoretically  and  experimentally.  As  is  shown  in  [3,  4],  in  this  case  the  viscosity  and  capillarity 
forces  are  not  considerable.  Therefore,  from  point  of  view  of  the  similarity  theory  the  dimensionless 
coefficient  of  the  gas-supply  rate  into  the  cavity  Q must  be  determined  by  the  dependence: 

Q=F{o,Fr). 


(11) 
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The  calculation  method  of  Q for  the  second  type  of  the  gas-leakage  was  proposed  for  the  first  time  by 
R.N.Cox  and  W.A.Clayden  [13]. 

Now  the  semi-empirical  formula  by  L.A.Epshtein  [3]  is  supposed  the  most  reliable: 


Here,  cx0  is  the  coefficient  of  the  cavitator  drag  when  o = 0: 


Cjr0  pv:sn 


where  Fx  is  the  cavitation  drag;  Sn  = n I)]  / 4 is  the  cavitator  area.  For  the  disk  cavitator  cx0  =0.82. 
Also,  the  more  exact  semi-empirical  formulae  were  obtained  taking  into  account  the  cavitator 
immersion  depth  H 

and  density  of  the  gas  supplied  into  the  cavity  pg , for  example  [6]: 


Q = 


(Fr-  1.35)o-‘ 75  Fr2(a3Fr4  -2.38c  J 25  ’ 


(14) 


where  H = H / Dn  . Computations  by  the  formula  (14)  give  satisfactory  coincidence  with  experimental 
data  when  Fr  = 4.8  -i-20  . 

The  experimental  dependencies  Q(g,  Fr)  showing  the  influence  of  Froude  number  on  the  gas 
leakage  by  vortex  tubes  are  presented  in  Fig.5  [3], 

The  gas-leakage  type  by  vortex  tubes  is  characteristic  also  for  the  case  of  supercavitation  flow  around 
hydrofoils  with  small  aspect  ratio. 

A criterion  of  transition  from  the  first  type  of  the  gas  leakage  to  second  type  is  given  by  the 
empirical  Campbell-Hilbome’s  criterion  [14]: 

<jFr  = l ; Fr  = 5 + 25  . (15) 


Values  of  a Fr  > 1 correspond  to  the  first  (portion)  type  of  the  gas-leakage,  values  G Fr  <\ 
correspond  to  the  second  type  of  the  leakage  (by  vortex  tubes).  As  the  experiments  show,  both  the 
mechanisms  of  the  gas  leakage  from  the  cavity  can  act  at  range  of  intermediate  values  of  G Fr  ~\. 
Such  transient  regimes  of  the  gas  leakage  are  very  unstable. 

A typical  experimental  dependence  of  the  gas  rate  coefficient  Q on  the  cavitation  number  G obtained 
at  testing  a series  of  cones  with  angles  from  30°  to  180°  (disk)  is  shown  in  Fig.  6 [3].  In  this  case  the 
value  of  the  Froude  number  was  Fr  = 13.3.  In  the  expressions  for  Fr  and  Q , so  called  "universal  linear 
dimension"  /J,  = Dn  jcx0  was  used  instead  Dn  to  analyze  the  experimental  data.  As  a result,  the 
relation  (11)  does  not  depend  on  cx0,  and  all  the  experimental  points  for  various  cones  lie  on  the  only 
curve. 


The  right  part  of  the  graph  1 corresponds  to  the  first  portion  type  of  the  gas  leakage,  the  left  part  of 
the  graph  2 corresponds  to  the  second  type  of  the  leakage  by  vortex  tubes,  the  middle  part  of  the 
graph  3 corresponds  to  the  transient  regime.  The  value  a = 0.081  corresponds  to  a boundary  between  the 
types  (15)  in  this  graph. 
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1.3.3  The  third  type  (gas-leakage  from  pulsating  supercavities) 


Denominator  in  the  formula  (13)  must  not  be  equal  to  zero.  It  follows  from  here  that  it  is  impossible  to 
obtain  the  cavitation  number  lower  than  some  minimal  value  at  any  increase  of  the  gas-supply  rate  into 
the  cavity: 

(16) 

In  practice,  when  the  gas  rates  are  very  great,  the  cavity  stability  is  lost  and  cardinal  change  of  the 
mechanism  of  the  gas  leakage  from  the  cavity  occurs  [8,  9], 


Theory  of  linear  stability  of  the  gas- filled  axially  symmetric  cavities  was  developed  by  E.  V.Paryshev 
[15].  Analysis  showed  that  the  dynamic  properties  of  gas-filled  supercavities  are  mainly  determined  by 
the  dimensionless  parameter 


P 


The  cavity  has  a series  of  fundamental  reduced  frequencies 


k=^  = 2nn 

V„ 


n = 1,2,..., 


(17) 


(18) 


where  ft)  is  the  circular  oscillation  frequency;  Lc  is  the  cavity  length.  The  following  value  of  the 
parameter  /)  corresponds  to  each  fundamental  frequency 


P„  - 1 + 


(nri)2 

6 


(19) 


Supercavities  are  stable  when  1 < /3  < /3,  = 2.645  and  unstable  when  /3  > 2.645. 

Thus,  the  similarity  parameter  /3  > 1,  that  is  equal  to  ratio  of  vapor  cavition  number  at  given  velocity 
and  its  real  value,  is  the  basic  similarity  parameter  for  cavitation  flows  together  with  the  cavitation 
number  (1)  and  the  Froude  number  (3).  The  value  /3  = 1 corresponds  to  the  natural  vapor  supercavitation. 
When  the  parameter  /3  increases,  significance  of  elasticity  of  the  gas  filling  the  artificial  cavity  increases 
as  well. 

The  separated  self-induced  oscillations  of  the  cavity  are  established  as  a result  of  the  stability  loss 
when  the  supply  is  permanent.  In  this  case  the  gas-supply  occurs  owing  to  periodic  separation  of  great 
portions  of  the  cavity  (air  pockets)  (see  Fig.  7). 

The  instability  phenomenon  and  arising  the  self-induced  oscillations  take  place  both  for  3-D  (in 
particular,  axially  symmetric)  and  for  2-D  ventilated  cavities.  In  the  work  [16],  a simple  "kinematic" 
theory  of  gas-leakage  from  pulsating  cavities  was  proposed.  It  gives  satisfactory  correspondence  to  the 
experiment. 
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1.3.4  Gas-leakage  when  the  cavity  closes  on  a body 


The  mechanism  of  gas  leakage  from  the  cavity  when  the  cavity  closes  on  a body  is  very  complicate  and 
depends  on  number  of  factors  such  as: 

a)  cavity  perturbations:  its  floating-up,  wave  deformations,  natural  distortion  of  the  free  boundary; 

b)  body  parameters:  shape  in  the  closure  place,  surface  roughness,  vibrations; 

c)  closure  conditions:  angle  of  free  boundary  inflow,  presence  of  the  liquid  and  gas  jets,  radial  velocity  at 
the  closure  place. 

According  to  the  known  theoretical  schemes  of  the  supercavity  closure  [17]  the  cavity  may  close  on  a 
solid  body  in  the  following  ways  (see  Fig.  8): 

a)  Riabouchinsky  scheme.  The  cavity  smoothly  closes  on  a solid  surface-closurer  located  in  the  cavity 
end.  In  the  closure  place,  the  body  diameter  is  Dh  <DC. 

b)  Joukowski-Roshko  scheme.  The  cavity  closes  on  a cylinder  with  diameter  equal  to  the  diameter  of  the 
biggest  cavity  section  Dc . 


c)  Gilbarg-Efros  scheme.  The  cavity  closes  on  a cylinder  with  diameter  Db  < Dc  with  forming  the  ring 
reentrant  jet. 

The  experience  of  work  with  the  artificial  cavity  shows  that  the  gas  leakage  may  be  compared  with 
the  rate  of  fluid  supplied  into  the  cavity  with  the  reentrant  jet.  Physically,  this  may  be  explained  that  the 
fluid  inflowing  into  the  cavity  aerates,  mixes  with  the  gas  and  rejects  from  the  cavity  as  a mixture  of 
gas  and  water. 

Using  the  momentum  conservation  equation,  the  cross  section  area  of  the  reentrant  jet  when  the 
cavity  closure  is  free  is  equal  [7]: 


c.v  nLK 

4 4 


(20) 


In  the  case  of  the  cavity  closure  on  the  body,  the  section  area  of  the  ring  reentrant  jet  remains  the  same. 
It  was  established  experimentally  that  when  the  cavity  closes  on  a circular  cylinder  with  diameter  Dh , 
the  gas  leakage  is  proportional  to  a ratio  of  section  areas  of  the  middle  cavity  part  nD2c  / 4 and  one  of  the 
cylinder  vD\  / 4 : 


( 

Q =k2V„Sj  1 
V 


D; 


Dh<D, 


(21) 


where  k2  is  some  empirical  coefficient.  When  Dh=  0 , we  have  free  closure  of  the  cavity  by  the  first 
type  of  the  gas-leakage.  When  Dh  = D , we  obtain  the  cavity  closure  by  the  Joukowski-Roshko 
scheme  with  the  zero  gas  leakage  (see  Fig.  8,  b). 

Theoretically,  we  can  also  obtain  the  zero  gas  leakage  when  the  cavity  closes  on  the  body  by  the 
Riabouchinsky  scheme  (see  Fig.  8,  a)  when  Db  <DC,  if  the  cavity  and  body  curvatures  are  coordinated 
at  the  closure  point  a=ah  (Fig.  9). 

Practically,  it  is  impossible  to  obtain  the  zero  gas  leakage  at  high  magnitudes  of  We  and  Re  numbers. 
However,  it  may  be  reduced  to  the  minimal  value  caused  by  presence  of  unremovable  small-scale 
perturbations.  Blowing  or  suction  of  fluid  is  applied  at  the  closure  place  to  realize  the  stable  cavity 
closure  [18].  For  example,  the  works  [18,  19]  propose  a scheme  of  stable  closure  of  the  cavity  on  the 
liquid  jet  outflowing  from  the  slot  in  the  closurer  of  the  elliptic  shape. 
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Fig.  8.  Schemes  of  the  supercavity  closure  on  a body:  a - Riabouchinsky;  b - Joukowski-Roshko;  c - 

Efros-Gilbarg 


We  have  discovered  experimentally  the  interesting  phenomenon  of  hysteresis  at  the  control  of  the 
ventilated  cavity  closing  on  the  body  by  regulating  the  gas-supply.  This  effect  increases  with  increasing 
the  body  surface  curvature. 

The  control  of  ventilated  cavity  is  possible,  if  single-valued  dependence  of  the  cavity  length  on  the  gas 
supply  rate  exists:  Lc  =LC{Q).  If  we  consider  neighborhood  of  the  point  of  the  cavity  closure  on 
rectilinear  generatrix  of  the  body  y - yh  (Fig.  10),  we  can  note  that  the  closure  angle  a will  be  increased 
with  increasing  the  cavity  dimensions: 

a(Ll)<a(L2)<a(L3),  I,<Z,2<I3 

due  to  increasing  both  the  ratio  Dc  / yb  and  the  relative  cavity  dimensions.  In  this  case  the  velocity  on  the 
cavity  boundary  changes  insignificantly  due  to  the  cavitation  number  varies  weakly.  Thus,  the  gas- 
leakage  value  Q depends  mainly  on  the  local  angle  of  the  cavity  closure  on  the  body: 

Q ~ Q(oc). 


Dc 


Fig.  9.  Scheme  of  the  fluent  cavity  closure  on  a body 


Dependence  of  the  gas-leakage  on  the  cavity  length  Qi  (Lc ) is  schematically  shown  by  dotted  line  in  Fig. 
10  when  the  cavity  closes  on  the  rectilinear  generatrix  of  the  body  y = yb  ■ 

Qualitatively,  appearance  of  hysteresis  behaviour  may  be  explained  by  locating  an  additional  triangle- 
shaped body  on  the  rectilinear  generatrix  (it  is  dashed  in  Fig.  10).  The  cavity  length  increases  with 
increasing  the  gas-supply  from  the  value  . When  the  cavity  closes  on  the  frontal  surface  of  the  triangle 
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x,  + x2,  the  cavity  forms  bigger  angles  of  closure  than  at  closure  on  the  rectilinear  generatrix.  For 
example,  a(x2)  is  bigger  than  ones  when  the  cavity  closes  between  the  points  x,  and  x4.  We  have 
a(x2)  = a(x4)  only  in  the  point  x4 . Therefore,  the  gas  rate  Q2  f°r  maintenance  of  the  cavity  closing  in 
the  point  x2will  exceed  rates  when  the  cavity  closes  in  the  points  x2-s-x4  and  equal  to  the  rate  for 
maintenance  of  the  cavity  closing  in  the  point  x4 . Asa  result,  the  cavity  length  will  very  rapidly  increase 
from  Lc  = L2  to  LC=L4. 

With  decreasing  the  supply  from  the  value  Q2 , when  the  cavity  closure  point  is  on  the  back  side  of  the 
triangle  x2  x3 , the  closure  angle  will  be  smaller  than  one  when  the  cavity  closes  on  the  rectilinear 
generatrix.  When  the  cavity  closes  in  the  point  x2 , the  rate  Qx  will  be  smaller  than  one  when  it  closes 
between  the  points  x,  and  x2 . As  a result,  the  cavity  length  will  decrease  very  rapidly  from  the  value 
Lc  = L2  to  the  value  LC<LV 


Fig.  10.  Scheme  of  a hysteresis  when  the  cavity  closes  on  a body 


Thus,  a hysteresis  loop  appears  in  the  graph  Q{LC ) that  is  formed  when  the  gas  supply  into  the  cavity 
increases  or  decreases. 

The  hysteresis  property  of  the  function  of  the  gas  supply  into  the  ventilated  cavity  may  by  practically 
used  to  stabilize  the  cavity  closure  on  the  body  by  respective  designing  the  body  shape. 


2 Approximate  calculation  of  the  steady  axisymmetric  supercavity 


As  it  is  known,  main  properties  of  supercavitation  flows  are  well  described  by  the  potential  theory  of 
ideal  incompresssible  fluid.  The  respective  mathematical  problems  are  related  to  the  class  of  mixed 
boundary  value  problems  for  harmonic  functions.  The  flow  area  is  bounded  by  the  solid  and  free 
boundaries.  The  velocity  distribution  is  given  on  the  solid  boundaries  (cavitator  or  hydrofoil),  the 
constant  pressure  pc  is  given  on  the  unknown  free  boundaries.  For  stationary  problems,  the  velocity 
constancy  along  the  free  boundaries  follows  from  the  pressure  constancy: 

Vc  = K„ -Jt  +<3  = const . (22) 

In  this  case  the  free  boundaries  are  lines  of  flow. 
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Potential  models  of  developed  cavitation  flows  do  not  describe  a real  flow  in  the  cavity  closure  zone 
where  capillarity  and  viscosity  are  considerable.  According  to  the  known  Brilluene  paradox,  existence  of 
a closed  cavity  in  the  imponderable  flow  is  impossible  [1],  This  resulted  in  creation  of  a number  of 
schemes  of  "mathematical"  closure  of  the  supercavity:  by  Riabouchinsky,  Joukowski-Roshko,  Effos- 
Gilbarg  et  al  [17]. 


2.1  Formulation  of  problem 


We  consider  steady  axially  symmetric  flow  around  a supercavitating  disk  in  unbounded  imponderable 
fluid  with  the  velocity  at  infinity  V^.We  introduce  the  cylindrical  coordinate  system  Oxr  connected  with 
the  cavitator  (Fig.  1 1).  Let  r = R(x)  is  the  equation  of  the  cavity  generatrix. 


2 


Fig.  11.  Scheme  of  an  axisymmetrical  supercavity 


The  Laplace  equation  for  the  perturbed  velocity  potential  cp(x,r ) is  valid  for  all  the  flow  area: 


dx  dr  r dr 


(23) 


On  the  cavity  surface,  the  potential  <p  must  satisfy  both  the  boundary  condition  of  zero  normal  velocity 
(the  kinematic  boundary  condition) 


d(p 

dr 


d(p  \dR 
+—  — . 
. dxjdx 


r = R(x) 


(24) 


and  the  cavity  pressure  to  be  constant  p = pc  (the  dynamic  boundary  condition).  Applying  the  Bernoulli 
equation  on  the  liquid  contour,  we  can  write  the  dynamic  boundary  condition  in  the  form: 


' d(p 
dr 


+ 


d(p 

dx 


+ 2F 


d(p  _ 2A p 
dx  p 


r = R(x) , 


(25) 


where  A p-p^~pc-  The  third  condition  is  the  condition  of  the  perturbations  to  be  vanished  at  infinity 


(p  — > 0,  x ^ + r — ^ °°  (26) 

Owing  to  great  mathematical  difficulties  an  exact  solution  of  the  boundary  value  problem  (23)-(26)  has 
been  not  obtained,  in  contrast  to  two-dimensional  problems  on  supercavitation  flows,  it  is  impossible  to 
apply  the  power  mathematical  theory  of  analytical  functions  and  conformal  mapping  in  this  case.  Some 
progress  of  applying  the  theory  of  generalized  analytical  functions  to  axisymmetric  problems  was 
outlined  last  years  [20]. 

All  known  results  of  the  theory  of  axisymmetric  cavities  have  been  obtained  by  numerical  methods  or 
by  approximate  asymptotic  and  semi-empirical  methods  and  perturbation  method  as  well. 
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P.Garabedian's  works  [21]  are  very  important  and  imperishable  for  the  theory  of  axisymmetric  super- 
cavitation flows.  He  sought  a solution  of  the  axisymmetric  problem  for  stream- function  \j/{x,r)  and 
consider  it  as  the  plane  flow  perturbation.  As  a result,  he  obtained  the  asymptotic  solution  for  main 
dimensions  of  the  supercavity  past  a disk  when  a — » 0 : 


D 


c»=cx(0Xl  + cr)l  c,(0)=  0.827, 


(27) 


where  D„  is  the  disk  diameter;  Dc  is  the  biggest  cavity  section  diameter;  Lc  is  the  cavity  length. 


A comparison  of  the  asymptotic  formulae  (27)  with  both  the  experimental  data  and  the  "exact" 
numerical  calculations  showed  their  high  precision  and  physical  accuracy. 


2.2  Approximation  relations 


A combined  approach  using  reliable  physical  representations,  general  conservation  laws  (such  as  the 
momentum  and  energy  theorems)  and  generalized  experimental  data  proved  to  be  the  historically  first 
and  most  productive  in  practical  calculations.  G.V.  Logvi-novich  was  a founder  of  this  direction  of  the 
super-cavitation  flow  investigation. 

We  present  firstly  some  useful  relations  for  the  axisymmetric  supercavities  obtained  by  approximating 
both  the  experimental  data  and  the  "exact"  numerical  calculations 


For  main  dimensions  of  the  axisymmetric  cavity  (the  mid-section  diameter  Dc  and  the  cavity  length  Lc ) 
in  the  case  of  the  cavitation  number  is  small  and  the  gravity  influence  is  negligible,  it  was  obtained  semi- 
empirical  relations  of  the  same  structure  as  in  (27)  [2,  3]: 


D„ 


4 


(28) 


where  k = 0.9  + 1 .0,  A are  empirical  constants.  When  a is  not  very  small,  it  is  usually  accepted  A ~ 2. 

Formulae  (28)  are  valid  not  only  for  the  disk  but  for  any  blunted  cavitators.  For  a non-disk  cavitator 
Dn  is  a diameter  of  the  cavity  separation  line.  For  the  drag  coefficient  of  the  blunted  cavitators  when  a 
is  sufficient  small,  the  approximate  Reichardt  fonnula  is  valid  [1-4]: 

Cx(<7)  = Cx0(l  + <T),  0<<7  <1.2,  (29) 


where  cx0  is  the  cavitator  drag  coefficient  when  a = 0.  Value  cx0  ~ 0.82  was  established  experimentally 
for  the  disk  cavitator.  For  cavitators  in  the  form  of  a blunted  cone  with  cone  angle  120  °<  /3  < 180  °,  it  is 
possible  to  use  the  formula: 

cx0  =0.13  + V0.0036/3-0.1719,  (30) 

that  was  obtained  by  approximating  the  experimental  data  for  the  disk  and  cones  [3]. 

Analysis  of  photographs  of  supercavities  past  the  blunted  cavitators  when  a is  small  permits  to 
distinguish  three  characteristic  parts  of  the  cavity  (Fig.  11).  The  frontal  part  1 adjoining  to  the  cavitator 
has  length  of  order  of  Dn . It  is  characterized  by  strong  curvature  of  the  free  boundary.  The  next  main 
part  II  occupies  approximately  V*  of  the  cavity  length  Lc . Finally,  the  tail  part  lit  has  unstable  unsteady 
boundaries  at  the  moderate  magnitudes  of  o . It  is  impossible  to  determine  them  using  the  potential 
theory. 
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The  frontal  part  boundaries  have  great  curvature  and  are  determined  by  the  cavitator  shape.  They  do 
not  depend  on  the  cavitation  number.  For  the  disk  cavitator,  the  empirical  "1/3  law"  is  usually  used  [2]: 


R_ 

IT 


1 + 


3x 

T 


V 


" J 


R. 


■<3^-5, 


where  R(x)  is  the  current  cavity  radius,  Rn  =Dn/2  is  the  disk  radius. 


(31) 


The  1HM  of  Ukrainian  NAS  accumulated  great  experimental  material  on  the  free  motion  of  super- 
cavitating  models  with  the  disk  cavitators  in  wide  range  of  velocities  V = 50  4-1400  m/sec.  The 
corresponding  cavitation  numbers  are  a = 1.4  • 10 4 4-  10 2 . Natural  vapor  cavities  are  formed  for  all  the 
pointed  velocity  range.  They  are  axially  symmetric  along  the  whole  their  length,  i.e.  are  not  deformed 
under  the  gravity  effect. 

For  the  main  supercavity  part  past  the  disk,  we  obtained  the  semi-empirical  formula  [22]: 

R2  =3.659  + 0.847(3t-2.0)-0.236o-(3c-2.0)2,  3c  > 2.0  (32) 


where  R=R/R„,  x = x/R„.  In  this  case,  the  experimental  data  was  used  which  were  obtained  for 
natural  vapor  supercavities  in  range  of  the  cavitation  number  a = 0.012  4-  0.057.  Results  of  measuring  the 
profile  of  the  steady  vapor  cavities  corresponding  to  various  values  of  a are  presented  in  Fig.  12.  Fig. 
13  presents  a photograph  of  the  vapor  supercavity  closing  on  the  jet  out-flowing  from  the  model  nozzle. 
The  frontal  cavity  part  at  3c  < 2.0  is  described  by  formula  (31). 


Fig.  12.  Universal  contour  of  the  axisymmetrical  supercavity 


Fig.  13.  View  of  the  natural  vapor  supercavity  past  the  disk:  Dn=  15  mm,  V=  105  m/sec,  a = 0.0204. 


From  the  formula  (32)  one  has  the  following  expressions  for  both  the  cavity  mid-section  radius  and  the 
cavity  length: 


*>£.-  3.659  + ^, 

R„  V <7 


(33) 


y 4 ..  3.595 

4 =-£-  = 4.0  + 

4,  o- 


(34) 


The  functions  Rc(a ) and  Lc(a)  which  are  calculated  by  formulae  (33),  (34)  in  range  of  cavitation 
number  a=  0.01  4-  0.06  are  shown  by  solid  lines  in  Fig.  14.  The  experimental  data  of  Fig.  12  also  are 
plotted  there.  For  comparison,  the  same  functions  given  by  asymptotic  formulae  by  P.Garabedian  (27)  are 
plotted  by  dotted  lines. 
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a b 

Fig.  14.  Approximate  dependencies  Rc(o)  (a)  and  Lc(<j)  (b)  calculated  by  formulae  (33)  and  (34) 


An  advantage  of  approximation  formula  (32)  is  that  it  gives  explicit  dependence  of  the  cavity  shape 
on  the  cavitation  number. 


We  also  present  other  approximation  formulae  for  the  main  dimensions  of  stationary  axisymmetric 
super-cavities  past  cones  and  the  cavitation  drag  coefficient  of  cones.  They  were  obtained  based  on  the 
numerical  calculations  in  "exact  formulation"  [23,  24]: 
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where  Dn  is  the  cone  base  diameter. 

cx  = cx0  + (0.524  + 0.672a)<7, 


cr  In—, 

1 (7 


0<(7<0.25,  — <a<  — . 

12  2 


(35) 


(36) 


cx0  =0.5  + 1 .81(ai  - 0.25)  - 2 (a  - 0.25)2, 


— <a <— , 
12  2 


cx0  = a(0.915  + 9.5a). 


0<a<- 


where  an  is  the  cone  half-angle.  When  a -a  0.  the  formulae  (35),  (36)  transform  into  the  asymptotic 
relations  by  P.Garabedian  (27). 


2.3  Energy  approach 


According  to  G.V.Logvinovich  work  [2]  we  use  the  theorem  of  momentum  to  obtain  the  approximate 
equation  for  the  main  supercavity  part.  We  display  the  check  volume  of  liquid  (Fig.  15).  It  has  such 
boundaries: 

5,  is  a plane  perpendicular  to  the  symmetry  axis  and  located  far  up-stream  ; 

52  is  a plane  perpendicular  to  the  symmetry  axis  and  passed  through  an  arbitrary  point  of  the  cavity 
axis  before  the  biggest  cavity  section; 

53  is  a cylindrical  surface  with  the  axis  coinciding  with  the  cavity  symmetry  axis  and  with  radius 
R2  — ^ °° . Let  R(x)  is  the  radius  of  the  cavity  section  by  the  plane  S2,  so  that  S = nR2  is  the  cross- 
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section  area.  We  write  the  theorem  of  momentum  for  the  displayed  fluid  volume  in  projection  on  the  x- 
axis: 

J PiYoo  + vx)2  ds  ~ J pvl  ds  = J pZ,  ds  - j pels  - Spc  - Fx,  (37) 

s2-s  st  s,  s2-s 

where  Fx  is  the  cavitation  drag.  Integrals  with  respect  to  the  cylindrical  surface  S3  vanish,  since  the 
perturbed  velocity  reduces  as  1 /r ! when  r—>  °°  (like  a flow  of  a dipole). 


s3 


Fig.  15.  Scheme  of  application  of  the  theorem  of  momentum 


From  the  Bernoulli  equation,  we  obtain  the  pressure: 

p = p~-p  Kyx  ~^pv2r,  (38) 

and  use  condition  of  conservation  of  the  mass  of  fluid  flowing  through  the  sections  S,  and  S2  - S: 

\pv:ds=  jp(V„+vx)2ds. 

S | S2-S 

It  follows  from  this  that 


j vxds  = V^S. 

s,-s 


(39) 


Substituting  (38)  and  (39)  in  the  Eq.  (37),  we  obtain  the  relation  for  cavitation  drag: 


px  = 4 ps~  j 


P»x 


dx  + 


pvr 


ds , 


(40) 


s,-s 


s,-s 


where  Ap-  px-  pc. 


If,  in  particular,  we  draw  the  plane  S2  through  the  cavity  mid-section,  then  S = Sc  and  vr  = 0.  Flence, 
the  second  integral  in  the  Eq.  (40)  is  strictly  equal  to  zero.  It  is  possible  to  show  that  the  first  integral  in 
the  Eq.  (40)  is  small  compared  to  the  first  term  when  a — > 0 [2].  That  is  why,  we  also  cast  it  out 
introducing  a correcting  multiplier  k.  After  that,  the  law  of  the  supercavitation  drag  may  be  written  in  the 
form: 


Fx  = kApSc  = kaSc 


2 


(41) 
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Introducing  the  cavitation  drag  coefficient  of  the  cavitator  (13)  in  the  Eq.  (41),  we  obtain  the  formula 
(28)  for  the  supercavity  mid-section  radius: 


Experimental  check  and  comparison  with  numerical  calculations  showed  that  the  correcting  multiplier  k 
weakly  depends  on  a and  has  magnitude  in  range  0.9  4-  1.0.  It  is  usually  accepted  k = 1 for  practical 
calculations. 


Then,  we  transform  the  Eq.  (40)  to  the  differential  equation  of  the  supercavity  section  expansion. 
Neglecting  the  first  integral  in  the  Eq.  (40),  we  may  consider  the  flow  in  the  plane  S2  as  pure  radial  one. 
We  have  the  boundary  condition  for  the  mean  part  of  the  slender  supercavity  surface 

vx  - VJi'(x).  (42) 

The  radial  velocity  of  fluid  at  the  arbitrary  radius  r is  determined  from  the  condition  of  the  mass 
conservation: 


2 nRvr{R)  = 2nrvr(r), 

then  we  obtain  taking  into  account  the  Eq.  (42): 

Mr)  = V„—  ■ (43) 

r 

The  integral  in  the  Eq.  (40)  represents  the  kinetic  energy  of  the  flow  in  the  fluid  layer  Ek . It  is  equal  to 


2 

Ek=  f ^-ds  = npVJRR')2\nr 


r=R 


(44) 


One  can  see  that  this  integral  value  is  equal  to  infinity,  i.e.  infinite  energy  is  necessary  for  the  radial 
expansion  of  the  cavity  in  the  two-dimensional  problem.  That  is  why  we  accept  two  non-strict  but 
plausible  assumptions: 

1)  Let,  according  to  [1],  the  outer  radius  of  the  flow  area  R2  has  high  but  finite  magnitude.  In  other 
words,  we  suppose  that  all  perturbations  of  the  fluid  motion  are  concentrated  in  a ring  with  the  inner 
radius  R and  the  outer  radius  R2  ■ This  supposition  is  plausible,  because  the  velocity  perturbation  decay 

very  rapidly  (as  1 / r 3 ) in  the  three-dimensional  problem. 

2)  Let  the  logarithm  magnitude 

^ht^  (45) 


does  not  depend  on  the  coordinate  x.  This  assumption  is  plausible  for  the  mean  part  of  the  slender  cavity 
when  cx  — x 0. 


Then  the  relation  (40)  may  be  rewritten  in  the  form 


Fx  = ApS  + ^-pVj(S')2.  (46) 

4n 

Differentiating  the  Eq.  (46)  with  respect  to  S,  we  obtain  another  form  of  the  equation  of  the  supercavity 
section  expansion  that  does  not  contain  the  cavitation  drag: 


d2S  2nAp 
u — T = f-  ■ 

dx2  pV2 


(47) 
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Passing  to  the  fixed  coordinate  system  in  the  Eq.  (47)  with  the  independent  variable  substitution  x = Vjt , 
we  rewrite  it  in  the  form 

3JS(T,t)  __  2gAp  (48) 

3r  P 

The  Eq.  (48)  assumes  the  obvious  mechanical  explanation: 

The  cavitator  passes  the  fixed  "observation  plane"  at  the  time  t = x and  makes  an  orifice  with 
diameter  Dn  in  that  plane.  Considerable  radial  velocities  arise  on  that  orifice  boundaries.  Then,  that 
orifice  expands  on  inertia  overcoming  the  pressure  difference  Ap  up  to  its  maximal  diameter  Dc  at  the 
time  t = T + tc . After  that,  the  orifice  decreases  and  closes  again  at  the  time  t = x + 2tc . 

The  time  interval  tc  = Lc  / 2 if  is  called  "the  time  of  the  complete  cavity  expansion".  After  the  cavity 
collapse,  the  wake  remains  past  it  where  the  mechanical  energy  transmitted  by  the  cavitator  to  the  fluid 
disperses  into  heat. 

The  Eq.  (48)  allows  to  calculate  an  expansion  of  each  cavity  section  separately,  i.e.  it  is  an  expression 
of  the  independence  principle  for  expansion  of  the  slender  axisymmetric  cavity  sections  by 
G.V.Logvinovich  [2], 

A full  cavity  profile  at  fixed  t may  be  calculated  when  t-2tc  <r<t . In  this  case  x-coordinates  of  the 
cavity  sections  are  x(z)  = Vjt , and  the  cavity  length  is  equal  to  Lc  - 2 tcV^ . 

The  more  strict  analysis  of  the  Eq.  (40)  [2]  permits  to  transform  it  to  the  following  form  (in  the  fixed 
coordinate  system): 

K"  c)S 

-0-tc)—  + (Sr-S)  = 0,  (49) 

2 at 

where  k is  some  correcting  coefficient.  The  integral  of  the  Eq.  (49)  is 

5(t)  = 5c+C(t-tc)2/K. 

Determining  the  arbitrary  constant  C from  the  initial  condition  S(0)  = Sn , we  finally  obtain  the  law  of  the 
cavity  section  expansion: 


r,  s„) 

(,  0 

2/a 

1- 

i — ,j- 

i — 

{ s<  j 

k G j 

(50) 


Passing  in  the  Eq.  (50)  to  the  coordinate  system  related  with  the  cavitator  x = Vj  , we  obtain: 


( R1  It 

R(x)  = R ,T-  1 — *- 


R- 


ix 


(51) 


For  non-disk  cavitators,  the  coordinate  x = 0 corresponds  to  the  plane  of  the  cavity  separation. 

A comparison  with  experiments  showed  that  value  of  the  correcting  coefficient  is  K =0.85.  However, 
the  exact  enough  results  are  obtained  when  K = 1. 

ft  is  easy  to  see  that  the  shape  of  the  main  part  of  the  axisymmetric  supercavity  represents  an  ellipsoid 
of  revolution  when  K = 1,  and  one  is  very  close  to  it  when  K = 0.85. 
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Experience  of  calculation  shows  that  the  expressions  (31),  (51)  must  be  "matched"  on  the  boundary 
between  the  frontal  and  the  main  cavity  parts  when  x = xl  ~ Dn  . As  a result,  we  obtain  so  called  formula 
of  the  G.V.Logvinovich's  composed  cavity: 


R(x)=Rn  | 1+— | - 


x < x, . 


R(x)  = R, 


1 - 

\ _«l_Y 

, 2 he 


(52) 


li  J 


where  R , =R(xl)  is  the  "agreement  section"  radius.  When  x,  =Dn,  we  just  obtain  from  the  Eq.  (31) 
that  i?,=  1 .913. 

The  "universal  contour"  of  the  ellipsoidal  cavity  is  plotted  by  solid  line  in  the  Fig.  12.  It  was 
calculated  by  the  formula  (52)  when  k=  1.  A comparison  with  the  experimental  data  shows  good 
coincidence  along  the  whole  cavity  length  right  up  to  the  point  of  the  cavity  closure  on  the  body. 


2.4  Asymptotic  approach 


Great  aspect  ratio  of  supercavities  when  a is  small  gives  background  to  apply  the  theory  of  flow  around 
slender  bodies  [25,  26]  to  calculate  them.  Basing  on  this  theory,  the  asymptotic  theory  of  a slender 
axisymmetric  cavity  was  developed  [27-30],  The  basic  mathematical  tool  of  this  theory  is  asymptotic 
method  of  singular  perturbations  [26],  A solution  of  the  axisymmetric  boundary  value  problem  is  sought 
in  the  form  of  decomposition  by  the  small  parameter  e . Its  value  is  inverse  to  the  cavity  aspect  ratio. 

It  is  well  known  from  the  theory  of  slender  bodies  that  the  potential  decomposition  near  the  x-axis  is 
[26]: 

V.  dR2  , 

q>  = In  r + ... 

2 dx 

To  evaluate  the  order  of  smallness  of  terms  in  the  Eqs.  (23)-(25),  we  substitute  there  orders  of 
smallness  for  all  the  variables: 


r ~ £,  R~e,  x ~ 1,  (p~e2\ne,  A p~e. 

As  a result,  neglecting  the  terms  of  high  order  of  smallness,  we  come  to  the  "internal"  problem  of  the 
first  approximation: 


d-<p  +l  d(p _0 
dr2  r dr 

(53) 

r «hen 
dr  dx 

r = R(x), 

(54) 

dx  p 

r = R(x). 

(55) 

We  note  that  the  condition  (26)  proves  to  be  lost  in  the  internal  problem. 

We  verify  by  direct  substitution  that  the  general  solution  of  the  Eq.  (53)  satisfying  to  the  kinematic 
condition  (54)  is: 

q>(x)  = ^— In  r + C(x),  (56) 

2 dx 


where  C(x)  is  some  arbitrary  function  that  may  be  determined  by  matching  with  the  "external"  solution 
satisfying  to  the  condition  (26)  [26]. 
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Substituting  the  solution  (56)  to  the  dynamic  boundary  condition  (55),  we  obtain  the  equation: 

2 C'(x)  2Ap 


d2R2 


In  R + - 


a2  • 


dx2  V„  pVj 

We  introduce  a new  arbitrary  function  'V(x)  instead  C(x)  using  the  equality: 


(57) 


C'(x)  = - 


F d2R2 


2 dx2 

Then  the  Eq.  (57)  becomes  the  same  form  as  the  Eq.  (47): 

R(x)  2A p 


ln'f'(x) . 


d2R2 

dx2 


ln- 


V(x)  pV2 


(58) 


if  we  consider  that 


In 


V(x) 

R(x) 


- p = const. 


Thus,  both  the  energy  theory  and  the  asymptotic  theory  of  slender  axisymmetric  cavities  yield  the  same 
equation  in  the  first  approximation. 

The  more  strict  analysis  of  the  boundary  value  problem  (53)-(55)  permits  to  obtain  the  differential 
equation 


d^R2 

dx~ 


ln-  + |— | = 


d_R_ 

dx 


2Ap 

~pvl 


(59) 


having  the  higher  precision  e2ln£  than  the  Eq.  (58). 

In  the  work  [30],  its  solution  was  constructed  by  the  method  of  matched  asymptotic  decomposition. 

However,  the  equation  of  the  first  approximation  (48)  is  the  more  convenient  for  practical 
calculations.  It  must  be  integrated  when  t > x with  the  starting  conditions: 


S(T,T)  = : 


kD: 


dS(  t,t) 
dt 


— S o > 


(60) 


where  S0  is  the  starting  velocity  of  the  cavity  section  expansion.  It  is  determined  by  the  cavitator  shape, 
i.e.  by  Dn  and  cx0 , and  also  by  the  model  velocity  V and  does  not  depend  on  the  cavitation  number. 

The  Eq.  (48)  reminds  the  Newton  equation  of  motion.  The  constant  p may  have  certain  physical 
sense  of  the  additional  mass  coefficient  of  the  expanding  cavity  section  (so  called  the  ring  model  of  the 
slender  axisymmetric  cavity  [1,  30]).  This  assumption  allows  to  determine  the  starting  velocity  of  the 
cavity  section  expansion  [29]: 

dR2 


- = 2 


t = T. 


(61) 


dt  y knpp  ’ 

However,  for  practical  calculations  it  is  rational  to  determine  both  the  constant  p and  the  initial  velocity 
of  the  section  expansion  S0  in  such  way  that  the  semi-empirical  relations  (28)  are  fulfilled.  Integrating 
twice  the  Eq.  (48),  we  obtain: 


dt  p 

S = S( T,  T)  + SQ  (t  - T)  - ^ - (t  — T)2  , 
2 P 


(62) 


(63) 
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where  kx  = 2k  / p . Substituting  the  expressions  (62)  and  (63)  into  obvious  relations: 


S(r,T  + tc)  = 


Kir 


dS(T,T  + tc) 
dt 


= 0, 


and  taking  the  Eq.  (28)  into  account,  we  finally  obtain: 


So  = 


(64) 


The  Eq.  (48)  is  an  equation  of  the  ellipsoid  of  revolution,  i.e.  it  correctly  describes  the  shape  of  the 
main  cavity  part.  To  take  into  account  the  frontal  part  (31),  it  is  necessary  to  use  the  condition 

S(t,t)  -=kD\ /4  (here,  £>,  is  the  "agreement  section”  diameter)  instead  the  first  of  starting  conditions 


(60). 

Fig.  16  gives  a comparison  of  the  cavity  shape  calculated  by  formulae  (52),  (63)  and  (32)  when 
<7  = 0.02. 


If  (7  < 0.01,  all  three  curves  coincide.  The  cavity  contour  approaches  to  the  ellipsoidal  cavity  shape  (52) 
when  o decreases  and  practically  coincides  with  the  latter  when  a<  0.01. 


Fig.  16.  Supercavity  shape  calculated  by  formulae:  1- (32);  2 - (52);  3 - (63) 


3 Calculation  of  unsteady  supercavities 

3.1  Independence  principle  of  cavity  expansion  by  G.V.  Logvinovich 


In  50s  G.V.Logvinovich  proposed  and  grounded  so  named  principle  of  independence  of  the  slender 
axisymmetric  cavity  section  expansion.  It  may  be  treated  as  a partial  case  of  the  general  mechanical 
principle  of  plane  sections.  That  principle  consists  in  that  the  law  of  the  cavity  section  expansion  almost 
does  not  depend  on  the  previous  and  the  next  motion  of  the  cavitator.  That  law  is  determined  by  the 
instantaneous  velocity  of  the  cavitator  at  passing  through  this  section  plane  and  by  the  pressure  difference 
Ap  = /;tx  - pc  as  well.  The  independence  principle  of  the  cavity  expansion  is  of  very  importance  to 
research  the  unsteady  supercavitation  flows. 


3.2  Equation  of  an  unsteady  supercavity 


In  the  case  of  unsteady  axisymmetric  cavities  the  independence  principle  consists  in  that  the  law  of 
arbitrary  section  expansion  about  the  trajectory  of  the  cavitator  center  is  determined  by  the  cavitator 
instantaneous  velocity  at  passing  this  section  V(t),  and  also  the  pressure  difference  Ap  = - pc  that  can 
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depend  on  time  at  fluctuating  the  external  pressure  px{t)  or  the  cavity  pressure  pc(t)  as  well.  The 
mathematical  expression  of  the  independence  principle  is  the  equation  of  expansion  of  the  unsteady 
cavity  section  [29]: 


d2S(r,t)_  k\Ap{t,t) 
dt 2 ~ P 


x(t)-Lc(t)<^<x(t), 


(65) 


where  Ap(T,t)  = p„(£)  + Pi(t)~  pc(t)  ■ Here,  r < I is  the  time  of  the  section/;  formation;  x(t)  is  the 
current  absolute  coordinate  of  the  cavitator  (see  Fig.  17);  p}(t)  is  the  perturbation  of  the  ambient 
pressure.  The  hydrostatic  pressure  px  can  change  from  a section  to  section  in  case  of  the  vertical  cavity 
in  ponderable  fluid. 


x-Lc(t) 


According  to  the  independence  principle,  the  equations  of  expansion  of  the  steady  cavity  (48)  and  the 
nonsteady  cavity  (65)  are  identical.  The  constant  /c,  = 2nl  p and  the  initial  velocity  of  the  section 

expansion  S(t,t)  = S0  have  the  same  magnitudes  (64).  In  this  case  V = V(r)  is  the  instantaneous 
cavitator  velocity  at  the  time  of  the  section  g formation. 

We  introduce  the  dimensionless  variables  by  formulae: 

x = Dnx\  Lc=DnL'c , S = D2S\  t = P = p V2  //. 


The  Eq.  (65)  for  the  dimensionless  variables  is  (the  primes  are  omitted): 

: , = -kiM?, t)=  — L-^ — 1 ^ (t),  t-Lc(t)< r < t. 

dt  2 

The  Eq.  (66)  must  be  integrated  with  dimensionless  initial  conditions: 

_ K dS  k ,A  / x / 

S=T-  TP—v{t)^-  ,=z- 


(66) 


(67) 


The  cavitation  number  is  a time  function,  if  the  external  pressure  p„(t)  and/or  the  cavity  pressure  p(  ( t ) 
fluctuate.  The  cavitator  velocity  V also  is  a time  function  in  the  general  case.  For  example,  in  case  of  the 
vertical  motion  of  the  supercaviating  model  under  action  of  only  cavitation  drag  the  cavitator  velocity 
and  the  cavitator  position  are  determined  by  the  equations: 


dV 

m = 

dt 


S„c 


xO 


1 + 


Patm  + PgH it)~  Pc(t) 
pV2/  2 


H{t)  = J V (s)ds, 
0 


where  m is  the  body  mass;  patm  is  the  atmospheric  pressure;  H is  the  cavitator  immersion  depth. 
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Twice  integrating  the  Eq.  (66)  with  taking  into  account  the  initial  conditions  (67),  we  obtain  the 
formula: 


S( r,  t)-  — + — — —\At]cx0  (t-r)  - 2V(t)  f(t-  u)o(it)du 

^ 4 T 

The  cavity  length  Lc(t)  is  determined  for  each  moment  by  the  equation: 


(68) 


S(t-Lc{t\t)=^.  (69) 

An  application  of  formula  (68)  has,  naturally,  restricted  limits  of  applicability.  However,  the 
experimental  tests  justify  its  adequacy  for  a wide  range  of  parameters. 


3.3  Experimental  confirmation  of  the  independence  principle 


The  experimental  test  of  the  independence  principle  was  carried  out  repeatedly  and  always  gave  good 
results.  So,  it  was  shown  in  the  specially  performed  experiments  [2]  that  in  case  of  the  model  impact 
acceleration  or  stoppage,  the  cavity  sections  remote  from  the  cavitator  expanded  some  time  in  the  same 
way  as  if  the  uniform  motion  continued. 

Examples  of  calculation  of  the  cavity  shape  at  the  variable  cavitator  velocity  and  also  along  the 
curvilinear  trajectory  are  given  there  as  well.  The  work  [30]  presents  results  of  experimental  testing  the 
calculation  of  the  unsteady  cavity  past  the  disk  with  periodically  varying  angle  of  slope.  A comparison  of 
the  experimental  and  calculated  cavity  shape  at  the  vertical  water  entry  of  bodies  is  given  there  as  well. 

A number  of  methodical  experiments  on  testing  the  independence  principle  was  carried  out  at  the 
IITM  UNAS  including  experiments  at  very  high  velocity  of  the  models  [31]. 

Fig.  1 8 shows  a sequence  of  motion-picture  frames  of  the  deformation  of  the  cavity  caused  by  rapid 
reducing  the  cavitator  drag  without  changing  the  cavitator  diameter.  The  experiment  was  performed  in 
the  hydro-tunnel  at  the  free-stream  velocity  I7,  = 8.9  m/sec,  the  diameter  Dn  of  the  cavitator  with 
variable  drag  [22]  was  equal  to  20  mm.  In  the  frames,  a well  expressed  "stepped"  boundary  separating 
the  initial  and  new-formed  cavity  parts  is  observed.  This  boundary  moved  with  velocity  in  full 
accordance  with  the  independence  principle. 

We  confirmed  the  independence  principle  validity  at  the  very  high  velocity  of  motion  in  the  special 
experiments  on  passage  of  the  supercavitating  model  through  a solid  steel  sheet  with  thickness  1.5  mm 
and  through  a hollow  obstacle  filled  by  air  as  well. 

Fig.  19  shows  a sequence  of  motion-picture  frames  of  the  process  of  the  model  passage  through  the 
thin  steel  sheet.  The  motion  velocity  is  V=  1000  m/sec,  the  cavitator  diameter  is  Dn  = 1.3  mm.  One  can 
see  in  the  frames  that  the  presence  of  the  obstacle  preventing  the  water  motion  in  longitudinal  direction 
did  not  influence  on  the  cavity  development  past  the  obstacle.  This  testifies  that: 

1 ) the  fluid  motion  near  the  slender  supercavity  occurs  predominantly  in  radial  direction; 

2)  the  adjoining  supercavity  sections  separating  by  the  steel  sheet  develop  in  such  way  as  in  the  case 
of  a flow  without  the  sheet; 

3)  the  longitudinal  gas  motion  together  with  the  supercavity  is  absent.  If  not,  it  would  cease  at 
passing  of  the  model  through  the  steel  sheet.  It  could  influence  on  the  supercavity  development  just  past 
the  obstacle. 

The  flow  directed  to  the  model  motion  direction  was  observed  only  after  the  cavity  closure. 
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Fig.  18.  Supercavity  shape  at  rapid  reducing  the  drag: 
through  the  obstacle: 

Dn  = 20  mm,  V„  =8.9  m/sec 


Fig.  19.  Passing  of  supercavitating  model 
Dn  = 1.3  mm,  V=  1000  m/sec 


The  experiment  on  passing  of  the  supercavitating  model  through  the  hollow  air  space  is  especially 
demonstrative  one.  The  obstacle  represented  a box  with  the  frontal  and  back  walls  made  of  steel  sheets 
having  thickness  1 mm  and  the  side  walls  made  of  glass  having  e thickness  10  mm  (Fig.  20).  In  the 
experiment,  we  used  the  model  with  the  length  280  mm  and  the  cavitator  diameter  Dn  = 3.3  mm. 

A sequence  of  frames  of  shooting  the  experiments  is  given  in  Fig.  21.  Passing  through  the  complex 
obstacle  containing  the  air  space  between  the  solid  walls  excluded  the  transition  of  both  the  pressure 
pulse  and  the  velocity  pulse  in  longitudinal  direction.  However,  the  cavity  developed  some  time  before 
and  after  the  hollow  box  as  if  it  was  absent. 

The  presented  sequences  of  frames  visually  demonstrate  validity  of  the  principle  of  independence  of 
the  cavity  section  expansion  and,  hence,  adequacy  of  the  approximate  calculation  method  based  on  that 
principle. 


4 Peculiarities  of  calculation  of  ventilated  cavities 


4.1  Similarity  parameters  in  the  artificial  supercavitation 


It  follows  from  the  above  that  the  basic  similarity  parameters  in  artificial  supercavitation  are  the 
cavitation  number  (1),  the  Froude  number  (3)  and  the  dynamic  parameter  /3  (17).  The  cavitation  number 

magnitudes  in  range  1 0 2<  <7  <0.1  correspond  to  this  type  of  flow. 

The  Froude  number  Fr  characterizes  the  distorting  effect  of  the  gravity  on  the  cavity  shape. 
Estimations  show  that  it  is  considerable  when  o Fr  < 2 [2], 
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Fig.  20.  Sketch  of  the  hollow  obstacle  Fig.  21.  Passing  of  supercavitating  model 

through  the  air  space:  Dn  = 3.3  mm,  V=  550  m/sec 


The  parameter  ft  > I characterizes  the  significance  of  elasticity  of  the  gas  filling  the  ventilated 
cavity.  It  is  the  basic  similarity  parameter  at  investigation  of  stability  and  unsteady  behavior  of  gas-filled 
cavities  [15].  When  the  parameter  (3  increases,  the  gas  elasticity  significance  increases  as  well.  The 

ventilated  cavities  are  stable  when  1<  /3  < 2.645  and  unstable  when  [3  > 2.645. 

In  practice,  it  is  impossible  to  attain  simultaneous  equality  of  all  the  numbers  <7 , Fr,  [3  for  the  full- 
scale  model  and  the  bench  testing  one.  That  is  why,  investigation  of  scale  effects  (i.e.  influences  of 
different  deviations  from  the  similarity)  is  of  importance  in  modeling  the  supercavitation  processes. 

We  enumerate  some  possible  causes  of  the  scale  effect  at  the  physical  modeling  of  the  artificial 
cavitation  flows. 

1)  Influence  of  the  flow  boundaries.  The  bench  testing  of  the  supercavitation  flows  usually  are 
carried  out  in  hydro-tunnels  and  hydro-channels  in  condition  of  the  bounded  How.  It  is  necessary  to  take 
into  account  the  influence  of  the  flow  boundaries  when  the  experimental  data  are  recalculated  to  the 
motion  in  unbounded  fluid. 

The  influence  of  the  flow  boundaries  on  the  supercavity  shape  and  dimensions  has  been  well  studied. 
The  solid  boundary  closeness  results  in  increasing  the  cavity  dimension,  and  the  free  boundary  closeness 
results  in  decreasing  them  at  the  same  value  of  the  cavitation  number. 

The  influence  of  the  flow  boundaries  on  the  gas  rate  coefficient  is  of  importance  in  the  case  of 
ventilated  cavities,  ft  was  shown  experimentally  that  the  approach  to  the  free  boundary  results  in 
decreasing  the  gas-leakage  rate  at  the  same  values  of  a and  Fr  [3,  10].  The  influence  of  the  cavity 
immersion  depth  on  the  gas-leakage  coefficient  at  the  second  type  of  the  gas-leakage  was  taken  into 
account  in  the  empirical  formula  (14). 

2)  Influence  of  the  gas  flow  in  the  cavity.  The  experience  shows  that  in  the  most  of  cases  the 
influence  of  gas  flow  within  the  artificial  cavity  may  be  neglected.  For  steady  cavities,  it  becomes 
considerable  only  in  that  cases  when  the  gas  flows  in  a narrow  clearance  between  the  body  and  the  cavity 
wall.  According  to  the  Bernoulli  equation  the  pressure  in  the  clearance  reduces,  and  the  cavity 
boundaries  may  be  deformed  negatively.  In  this  case  local  deformation  of  the  cavity  wall  may  be 
approximately  calculated  by  the  perturbation  method  [32]. 
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3)  Supercavities  in  a bubble  flow.  A peculiar  scale  effect  may  arise  in  the  case  of  the  gas-liquid 
bubble  flow  around  the  supercavitating  bodies.  We  showed  experimentally  and  theoretically  that  the 
supercavities  can  absorb  gas  from  the  two-phase  bubble  flow  [33]. 

When  the  gas  concentration  is  sufficient,  this  process  causes  to  the  considerable  gas  supply  into  the 
cavity  and,  as  a result,  to  increasing  the  cavity  dimension. 


4.2  Taking  into  account  of  gravity  effect 


As  was  said  above,  the  gravity  effect  on  the  horizontal  supercavity  consists  in: 

1)  the  cavity  axis  deformation  upwards  (floating-up  of  the  cavity  tail); 

2)  the  cavity  section  shape  deformation. 

The  cavity  axis  deformation  in  ponderable  fluid  is  the  main  and  greatest  perturbation  of  the  cavity. 


Fig.  22.  View  of  the  ventilated  supercavities:  = 8.9  m/sec,  Fr  = 24.5;  a - (7  = 0.0334,  b - <J  = 

0.0644 


Fig.  22  demonstrates  photographs  of  ventilated  cavities  at  the  same  Froude  number  and  the  different 
cavitation  numbers.  The  first  form  of  the  gas  leakage  from  the  cavity  was  observed  in  both  cases. 

In  the  works  [31,  33]  these  deformations  have  been  investigated  in  detail  by  the  perturbation  method. 


Fig.  23.  Scheme  of  a supercavity  in  ponderable  fluid 
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However,  for  practical  calculations  it  is  usually  enough  to  evaluate  the  axis  deformation  in  the  first 
approximation  by  the  theorem  of  momentum  [2],  The  calculation  scheme  is  presented  in  Fig.  23. 

We  consider  the  balance  of  the  cavity  length  element  with  the  radius  R(x).  Its  momentum 
pnR2V^hg  must  be  equal  to  the  buoyancy  force  pgQ . Thus,  we  have: 


K(x)=- 


gQ(x) 


nV„R\x) 

Integrating  this  equation  along  the  x-axis,  we  obtain 

gCO 

*2C0 


K F { 


ds . 


(70) 


Here,  Q(x)  is  the  cavity  part  volume  from  0 to  x,  R(x)  is  the  current  cavity  radius. 

In  the  work  [33],  the  simple  approximate  formula  was  obtained  for  the  axis  deformation  caused 
to  the  gravity  effect: 


k(x)-- 


+ 

3 Fr2 


Fr,= 


F 


(71) 


Here,  Frt  is  the  Froude  number  with  respect  to  the  cavity  length.  Formula  (71)  may  be  used  in  ranges 
0.05  < (7  <0.1,  2.0  <Fr,  <3.5.  Fig.  24  gives  a comparison  of  calculation  by  formulae  (70)  (curve  1) 

and  (71)  (curve  2)  with  the  experimental  data  [30],  Different  marks  correspond  to  different  experimental 
conditions. 


Fig.  24.  A supercavity  axis  deformation  in  ponderable  fluid 


Deformation  of  the  cavity  section  shape  in  the  cross  gravity  field  was  investigated  by  the  perturbation 
method  in  the  works  [30,  32].  The  calculations  showed  and  the  experiments  confirmed  that  a liquid  crest 
with  top  directed  up  into  the  cavity  is  formed  at  the  bottom  cavity  wall  in  this  case.  The  crest  height 
increases  from  the  cavitator  to  the  tail.  Sometimes,  this  crest  is  small  even  at  the  cavity  end.  In  other 
cases  the  crest  top  may  attain  the  upper  cavity  wall  and  divide  cavity  on  two  tubes  just  past  the  mid- 
section (see  Fig.  25) 

A criterion  of  the  gravity  effect  on  the  supercavity  was  obtained  by  theoretical  way  in  the  work  [32]: 

v=o4oFr2>  1.5.  (72) 

When  this  condition  is  fulfilled,  the  cavity  is  not  destroyed,  and  the  water  crest  does  not  close  up  with  its 
upper  wall  even  at  the  cavity  end.  With  approaching  the  parameter  v to  the  limiting  magnitude  V = 1.5, 
two  vortex  tubes  become  the  more  expressed  at  the  cavity  end. 
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A comparison  with  the  experiment  showed  that  the  criterion  (72)  is  the  more  exact  than  the  Campbell- 
Hilbome  criterion  (15).  Indeed,  we  have  a Fr  =0.82  < 1;  v = 3.73  > 1.5  for  the  experiment  in  Fig.  22,  a. 
For  the  experiment  in  Fig.  22,  b,  we  have  a Fr=  1.58  > 1;  V = 9.82  > 1.5.  Thus,  in  the  first  case  the 
Campbell-Flilbome  criterion  does  not  allow  confidently  to  determine  the  gas-leakage  type. 
Simultaneously,  the  criterion  (72)  determines  the  portion  gas-leakage. 

It  is  possible  to  determine  the  following  intervals  of  changing  the  parameter  v that  correspond  to 
different  levels  of  perturbation  caused  by  gravity: 

1 )  interval  V =1+2  is  characterized  by  high  level  of  perturbation,  when  the  cavity  past  the  mid- 
section transforms  into  the  vortex  tubes; 

2)  interval  V = 2 + 4 corresponds  to  the  moderate  or  considerable  level  of  perturbation.  In  this  case 
the  water  crest  height  is  lower  than  the  unperturbed  cavity  radius,  although  may  be  close  to  it. 

3)  interval  V =4+10  is  characterized  by  weak  level  of  perturbation.  The  gravity  effect  may  be 
neglected  at  all  when  v > 10. 


Fig.  25.  Deformation  of  a ventilated  supercavity  axis  in  ponderable  fluid 


Fig.  26  represents  the  calculation  by  the  perturbation  method  of  the  cross  cavity  section  shapes  when  a = 
0.06,  Fr  = 10.0  [33].  In  this  case  v = 1.47.  The  dimensionless  coordinate  x = 2x/Lc  was  used  at  the 
calculation.  In  the  sections  3c  = 0.25  and  x = 0.50  the  gravity  effect  does  not  yet  become  apparent,  in  the 
section  3c  = 0.75  the  small  compression  from  the  bottom  is  appreciable.  In  the  cavity  mid-section  3c  =1 .00 
the  section  shape  differs  from  the  circular  one.  When  x = 1.25,  the  deformation  is  considerable.  On  the 
bottom  it  is  close  to  the  half  of  the  cavity  radius.  The  crest  directed  to  up  is  well  visible.  The  cavity  is 
destroyed  in  the  next  sections. 


Fig.  26.  Gravity  effect  on  the  cross  cavity  section  shape:  a = 0.06,  Fr  = 10.0 
1 - 3c  = 0.25,  2 -3c  =0.50,  3 - 3c  = 0.75,  4-  3c  =1.00,  5-  3c  = 1.25 
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4.3  Taking  into  account  of  angle  of  attack 


The  transverse  force  F arising  on  the  cavitator  non-symmetric  about  the  free-stream  results  in  the 

cavity  axis  deformation  (Fig.  27).  According  to  the  theorem  of  momentum,  the  impulse  of  the  transverse 
force  on  the  cavitator  must  correspond  to  change  of  the  momentum  in  the  wake  that  is  the  same  by 
magnitude  and  opposite  by  direction.  This  means  that  if  the  force  on  the  cavitator  is  directed  to  up,  then 
the  cavity  axis  must  be  deflected  to  down  (see  Fig.  27)  and  vise  versa. 

Applying  the  theorem  of  momentum  allows  to  estimate  easily  this  bend  [1].  The  impulse  F t must  be 
equal  to  momentum  of  the  cavity  length  unit  - pnR2Vjhv . Hence,  we  have: 


hf(x)  = 


Fy 

pnV2R2(x) 


Fig.  27.  Scheme  of  a supercavity  past  an  inclined  cavitator 


Integrating  this  equation  along  the  x-axis,  we  obtain 

*/(*)  = 


npVl  Jo^2('?) 


A 

1 


ds 


(73) 


Components  of  the  force  acting  on  the  cavitator-disk  inclined  to  the  stream  with  the  angle  a can  be 
approximately  calculated  by  formulae  [1]: 


Fx  = Fx0  cos~  a,  F = Fx0  sin  a cos  a, 


(74) 


where  Fx0  is  the  disk  drag  when  a = 0.  The  fonnulae  (74)  well  coordinates  with  experiment  when 
a < 50°.  They  may  be  applied  not  only  for  a disk  but  also  for  blunted  cavitators  of  other  shapes. 


Fig.  28.  A superavity  axis  deformation  past  the  inclined  cavitator:  a =0.08,  a =33.5° 


An  analysis  of  the  cavitator  orientation  effect  on  the  cavity  shape  using  the  perturbation  method  was 
given  in  the  work  [34].  As  a result,  the  simple  approximate  formula  was  obtained  to  determine  the 
cavity  axis  deformation  under  action  of  the  lift  on  the  cavitator: 

pV2nDl  ’ 


hf(x)=-cyRn  (0.46-cr  + x\ 


(75) 
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where  x = 2x/Lc.  Fig.  28  gives  a comparison  of  the  calculation  by  the  formula  (75)  with  the 
experimental  data  [3 1], 

Comparing  the  formulae  (74)  and  (78),  we  can  see  that  the  cavity  axis  bend  caused  by  the  gravity 
increases  according  to  the  quadratic  law,  and  one  caused  by  the  lift  on  the  cavitator  increases  according 
to  the  linear  law.  Therefore,  it  is  impossible  to  neutralize  the  gravity  effect  using  the  inclined  cavitator 
along  the  whole  cavity  length. 

In  skewed  stream,  the  cross  cavity  sections  also  are  deformed.  Since  a projection  of  the  inclined  disk 
on  the  vertical  plane  is  an  ellipse,  then  the  cross  cavity  sections  before  the  mid-section  will  be  elliptic 
with  ratio  of  the  half-axes 


— = cosa. 

K 


(76) 


In  this  case,  the  average  radii  in  the  middle  cavity  part  before  the  mid-section  are  approximately  equal  to: 

R{x)=  R0(x)cosa , (77) 


where  R0  is  the  radius  of  axisymmetric  unperturbed  cavity  [2], 

The  simple  formulae  (71),  (75)  may  be  used  for  approximate  calculation  of  both  the  stationary 
and  non-stationary  supercavities  (using  the  independence  principle).  They  are  very  convenient  for  using 
in  codes  for  computer  simulation  of  the  cavitation  flows. 


4.4  Equation  of  the  mass  of  gas  in  the  cavity  balance 

The  mathematical  model  of  the  unsteady  ventilated  cavity  must  include  the  equation  of  the  mass  of  gas  in 
the  cavity  balance  [15]: 

dm 

— — = to,„  -mont,  (78) 

at 

where  mg(t)  = p (t)Q(t ) is  the  mass  of  gas  in  the  cavity;  (Q(t)  is  the  cavity  volume;  p„  is  the  density  of 
gas  in  the  cavity;  min  and  mout  are  the  mass  rates  of  the  gas  blowing  into  the  cavity  and  carrying  away 
from  the  cavity,  respectively.  Owing  to  the  water  thermal  capacity,  the  process  of  the  gas  expansion  in 
the  cavity  usually  may  be  considered  as  isothermal  one: 

pg(t)  = CpJt),  (79) 

where  C is  the  constant.  Taking  the  relation  (79)  into  account,  we  rewrite  the  equation  of  the  mass  of  gas 
in  the  cavity  balance  in  the  form: 

^-[Pc  (t)Q(t)]  = P„  [fi„  - Qout  col  (80) 

at 

where  Qin  is  the  volumetric  air-supply  rate  referred  to  the  ambient  pressure  p„;  Qoul(t ) is  the 
volumetric  air-leakage  rate  from  the  cavity.  The  cavity  pressure  pc{t)  is  assumed  to  change 
synchronously  along  the  cavity  length. 

Strictly  speaking,  the  Eq.  (80)  may  be  applied  when  the  parameters  change  not  very  rapidly.  Really, 
the  pressure  perturbations  pc (t)  extend  along  the  cavity  with  the  sound  speed  in  the  vapor-gas  medium 
filling  the  cavity  ag . Characteristic  time  of  changing  the  pressure  in  the  cavity  with  length  Lc  is  equal: 
Tj  = Lc  / a . Let  the  characteristic  frequency  of  the  unsteady  process  is  / Hz.  Then,  supposition  about  the 
synchronism  of  changing  the  pressure  along  the  cavity  is  equivalent  to  the  inequality  TJ  « 1 If.  From 
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here  we  derive  estimations  for  both  the  admissible  frequency  of  the  process  and  the  dimensionless 
Strouhal  number  (10): 

a„  a„ 

f«—,  Sh «—.  (81) 

L V 

C oo 

The  cavity  pressure  is  constant  in  the  case  of  the  natural  vapor  supercavity,  i.e.  a = const.  Therefore, 
the  Eq.  (80)  becomes  needless. 

In  practice,  we  usually  have  Qin  = const.  The  value  of  the  gas-leakage  from  the  cavity  Qout  depends 
on  both  the  cavitation  number  and  the  Froude  number  and  on  a number  of  other  factors  as  well  [3,  6].  It  is 
obvious  for  steady  ventilated  cavity  that  Qin  = Qout . 

The  function  Qin  (?)  is  given,  and  Q,ml  (?)  is  an  unknown  time  function  in  problems  of  the  ventilated 
cavities  control  by  regulating  the  gas  supply.  As  has  been  said  yet,  determining  the  form  of  function 
Qom  (0  f°r  the  concrete  flow  conditions  is  the  most  difficult  problem  of  the  artificial  cavitation  theory. 


Conclusion 

The  described  approximate  mathematical  model  of  the  unsteady  axisymmetric  cavity  is  very  convenient 
to  develop  fast  calculation  algorithms  and  codes  permitting  to  perform  computer  experiments  with 
dynamic  display  of  the  unsteady  cavity  shape  and  other  necessary  information  on  the  computer  screen 
immediately  during  run-time.  The  existing  methods  of  "exact"  numerical  calculation  of  the  unsteady 
supercavities  are  useless  for  these  purposes  due  to  their  awkwardness  and  long  run-time.  Accuracy  of 
the  obtained  results  is  established  by  direct  comparison  with  experiment. 

In  our  next  lecture,  results  of  computer  simulation  of  the  unsteady  supercavitation  processes  of 
different  types  will  be  presented. 
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